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Free-electron lasers (FELs) can now generate temporally short, high power x-ray pulses of un- 
precedented brightness, even though their longitudinal coherence is relatively poor. The longitu- 
dinal coherence can be potentially improved by employing narrow bandwidth x-ray crystal optics, 
in which case one must also understand how the crystal affects the field profile in time and space. 
We frame the dynamical theory of x-ray diffraction as a set of coupled waves in order to derive 
analytic expressions for the spatiotemporal response of Bragg scattering from temporally short in- 
cident pulses. We compute the profiles of both the reflected and forward scattered x-ray pulses, 
showing that the time delay of the wave r is linked to its transverse spatial shift Ax through the 
simple relationship Aa; = cr cot 6, where 9 is the grazing angle of incidence to the diffracting planes. 
Finally, we apply our findings to obtain an analytic description of Bragg forward scattering relevant 
to monochromatically seed hard x-ray FELs. 

PACS numbers: 41.50.+h, 41.60.Cr 



Free-electron lasers (FELs) based on self-amplified 
spontaneous emission (SASE) [3 [2] are a new source 
of x-rays whose brightness is many orders of magnitude 
larger than those of more traditional third-generation 
synchrotrons. Nevertheless, because SASE radiation is 
initialized (or seeded) by the shot noise of the electron 
beam, it comprises many temporal/longitudinal modes 
[3] . The bandwidth of this chaotic light is limited to that 
of the FEL process, which for high-gain FELs means that 
the normalized frequency bandwidth Aw/w ~ p, where 
p is the FEL Pierce parameter that is typically between 
10--^ and 10-3 for hard x-ray FELs. 

While SASE FELs are hmited to Auj/w ~ p, one can 
decrease the output bandwidth and increase the longi- 
tudinal coherence by seeding the FEL with coherent ra- 
diation from an external source. For hard x-rays where 
an external source is not available, it was proposed that 
one could improve the longitudinal coherence through 
self-seeding - by putting the SASE light through an x- 
ray monochromator one can obtain a narrow bandwidth 
source with which to generate longitudinally coherent 
FEL radiation. The primary difficulty with the scheme 
proposed in [Ij was that the monochromator delayed the 
radiation by several picoseconds, which in turn required 
a long 40 m) magnetic chicane to delay the electron 
beam by the same amount. One way of alleviating this 
is to use two separate electron bunches [S]. 

Recently, an alternative seeding method was proposed 
[51 [7] that takes advantage of the time response of Bragg 
forward scattering to generate the seed. In this scheme, 
dubbed the "wake monochromator," a single Bragg crys- 
tal is employed in transmission to generate a delayed, 
monochromatic signal (the wake) that is then used to 
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seed the FEL. Since the wake follows the SASE pulse by 
~10 fs, a modest chicane is sufficient to induce the re- 
quired delay in the electron beam. The physics discussed 
in [5] was in terms of a band-stop filter whose sharp spec- 
tral feature naturally gave rise to a long, monochromatic 
temporal signal that was delayed due to causality. Here, 
we give a more thorough analysis of the process, obtain- 
ing analytic expressions for the spatiotemporal shape and 
power of the monochromatic wake generated by an inci- 
dent SASE pulse. 



The temporal dependence of x-ray Bragg diffraction 
in crystals has been studied previously in several publi- 
cations 'H'-14' by Fourier transforming the frequency do- 
main Bragg diffraction amplitudes, which are well-known 
from the classical dynamical theory of x-ray diffraction 
[T5l419j . Here, our focus is on the spatiotemporal dy- 
namics that play an essential role in the efficacy of the 
FEL self-seeding. For this reason, the problem is solved 
directly in time and space using a system of space-time 
coupled wave equations. 



We begin Sec. |T]by deriving the coupled wave system 
relevant to symmetric Bragg scattering, and subsequently 
solve for the field profiles in both reflection and trans- 
mission. We show that the temporal profile consists of 
a sequence of power maxima whose location depends on 
the crystal extinction length A and grazing incidence an- 
gle 9 in reflection, while in transmission the time delay 
also varies inversely with the crystal thickness d. In ad- 
dition, the profiles are displaced along the transverse di- 
rection Xo by an amount proportional to the delay r: 
Axo = CT cot 6. Finally, we use some well-known FEL 
physics to apply the theory of forward scattering to the 
FEL self-seeding scheme, obtaining analytic results asso- 
ciated with the wake monochromator of Ref. [5] . 
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I. BRAGG SCATTERING IN THE TIME 
DOMAIN 

We are predominantly interested in understanding the 
basic physics and deriving a few simple analytic expres- 
sions for the time response of Bragg forward scatter- 
ing from temporally short incident and laterally confined 
pulses. To simplify the subsequent analysis and avoid 
notational complications, we assume that the crystal is 
symmetric, meaning that the crystal planes responsible 
for Bragg scattering are parallel to the crystal surface. 
This restriction is most relevant to the FEL monochro- 
mator we subsequently study. We take the optical axis 
of the incident radiation to make an angle 9 from the 
surface as shown in Fig. [l|a). Maxwell's wave equation 
for the electric field E in the crystal is 



1^ 



where J is the current density induced in the crystal 
by the radiation and c is the speed of light. The linear 
response of the medium is given by the crystal polariz- 
ability x(t). Denoting the polarization components by a 
superscript s, in the linear approximation for x-ray ener- 
gies far above any atomic resonances we have 



s' 
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where in the second line we have used the periodicity of 
the crystal to expand the polarizability x as a Fourier se- 
ries. We decompose the electric field vector as a sum of 
two orthogonal polarizations, E — E'^e.a+E'^e.-,,, with 
parallel and perpendicular to the scattering plane. In 
this case, the non-diagonal components (namely, Xh 
s' 7^ s) vanish and the polarization components decou- 
ple. For simplicity, we assume cr-polarization and subse- 
quently drop the polarization dependence. 

The electronic polarizability in ([2| strongly couples 
electromagnetic waves whose wave-vectors differ by /i, 
which results in strong Bragg reflection. We assume that 
this condition is satisfied for one reciprocal lattice vector, 
which in the symmetric geometry implies that h — —hz. 
Since the Xh ^-re assumed to be time independent, it is 
natural to take a Fourier transform with respect to i; 
however, our goal is to compute the temporal response, 
and we find it more convenient to remain in the time do- 
main. Instead, we introduce the slowly varying eikonal 
field amplitudes associated with two interacting waves, 
writing the electromagnetic wave as 
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(3) 



where cfep is the carrier frequency and (feg sin 9, 0, fcg cos 9) 
is the reference wave-vector. We assume that the field 
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FIG. 1: Geometry of Bragg scattering, (a) shows that the 
incident pulse travels along the optical axis that makes an 
angle 6 with respect to the crystal surface, (b) plots the 
geometry in the characteristic coordinate plane C'"C'i with 
the crystal entrance and rear surfaces located at = ^' (the 
line AQ) and = ^' — d (the line CP), respectively. We also 
include the integration contours ABQ and ACPQ that are 
used to obtain the solution for the reflected and transmitted 
waves. 



envelope functions have slow spatiotemporal variations 
with respect to the carrier frequency cfco, with 
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The reference carrier frequency cfcp can be freely chosen 
provided that Q hold, with any additional slow temporal 
dependence being accounted for by Eg^h] we choose ko to 
satisfy Bragg's condition: feg sin 9 = h/2. 

We have oriented our axes such that the scattering 
plane defined by the reciprocal lattice vector and the op- 
tical axis lie in the (x, z) plane, and we assume that the 
crystal is uniform in the {x,y) plane. Thus, the funda- 
mental field solutions are plane waves in the directions 
parallel to the crystal surface. For simplicity we will ne- 
glect the trivial y-dependence, with the general electric 
field E being written as the superposition 



Eo,h{x,z,t) ^ dq e"^''£o,h{q,z,t). 



(5) 



We now use the slowly-varying assumption Q to re- 
duce the second order wave equation ([T]) to two first-order 
equations for the field amplitudes Eq and Eh , an approach 
that was first developed in this context by Takagi pU[|21) . 
Since we assume that the crystal lattice is perfect, we find 
it convenient to use z and t as the relevant coordinates 
rather than the standard spatial coordinates along the 
wave-vectors (see, e.g., |18j il9j US]); thus, our ap- 
proach is similar to those of [Ml [25] , We insert the forms 
([5]) and ([3| into the wave equation (!]), and drop the term 
^ V(V • E) since the field remains approximately trans- 
verse in the crystal. By neglecting higher-order deriva- 
tives of the field amplitudes and matching fast phases. 
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we obtain the following set of slowly-varying equations: 
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where for convenience we define Xh = X-h- We have 
attached tildes to ao,/i because while they appear similar 
to the usual deviation from Bragg's condition ao.h, they 
have a slightly different definition and interpretation in 
the present context, ao represents the incidence wave- 
vector's difference from the vacuum condition due to g ^ 
0, while a/i is the deviation of the reflected carrier wave 
from Bragg's condition, with 
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our choice fco sin — h/2 implies that 



ao = a/i = — cos ( 
fco 



2fcQ 



(8) 
(9) 
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Note that the central wavenumber fco and the crystal pa- 
rameters xoj Xhi ^"^d Xh ^'^^ functions of angle 9. To 
solve we introduce the characteristic coordinates 



C = i (ct sin 9 — z) , ^ = ^ {ct sin 9 + z) , 

and the reduced field amplitudes Aq and Ah via 

So = e"^*-^«(""Xo/2)(C+«)/sine_^^ 

£h = e^^''"^"^^°/'^'>^'^+^^/ ''"''^ Ah- 
Then, the coupled wave system (|6])-([7]) reduces to 
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(12) 
(13) 
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Here, we see that if the coupling Xh ~^ Oj the forward- 
going wave Aq is a function of ( only, while the reflected 
wave Ah is a flxed function of i^. In the crystal, the two 
waves interact as shown, with each obeying the associ- 
ated second-order equation 
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where the extinction length A = (27r/fco) smO / ^/XhXh 
approximately independent of the angle 9 for any given 
reflection [this can be seen by using Xh ~ ^(X/i) ^ sin^ 9 
and fco = /i/(2sin6')]. 



The linear system (14) can be solved for specifled 
boundary conditions using the Riemann method as was 
done in |22) . and which is reviewed in, e.g., |19j and in 
the Appendix. For Bragg scattering, the reduced fields 



satisfy (14) subject to the constraints that the forward- 
going wave Ao is a prescribed function along the front 
crystal surface at z = while the reflected wave vanishes 
along the rear surface defined by z = rf. In terms of the 
characteristic coordinates shown in Fig. [ijb), we have 



(16) 



In the following two sections we solve the system ( 14 ) 
subject to ( |16[ ) assuming a Gaussian incident pulse that 
is temporally short and confined laterally in space. We 
will find relatively simple, approximate analytic expres- 
sions for the electric field when we can neglect multiple 
scattering of the waves off the crystal surfaces, meaning 
that we only consider the time interval following the main 
pulse that is smaller than {2d/c)/ sn\9. For a crystal of 
thickness d ~ 0.1 mm this time interval is ~ (300/ sin 0) 
fs. The field at longer times can be built up by inserting 
these solutions into the Riemann integrals in an itera- 
tive manner [551 IMI j but this is beyond the scope of the 
present study. 



A. Bragg diffraction: the reflected wave 

In the Appendix we use Riemann's method to show 
that on the front crystal surface the reflected fleld is given 
by 
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where ^inc is the initially incident wave and the Riemann 
function Rh for the reflected wave satisfles the adjoint 
equation associated with the (formally self-adjoint) sys- 
tem (fTil) 
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along with the auxiliary conditions 
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The Riemann function satisfying (18)-(19) is given by 
(A9) [H]; for the solution (17) we need Rh along the 



front crystal surface defined by the line AQ; here, we 
have ^' = with point Q located at ^ = C, for which 



RhiCX;C,C)^2 



Ji[27r(C-C0/A] 
27r(C-C')/A ■ 



(20) 



The physical electric field amplitude of the reflected 
wave can be obtained by applying the field definitions 
(fT2|)-(ll3|) to the solution (17) once we specify the inci- 



dent field. We will be interested in the response from 
a temporally short incident wave directed along the op- 
tical axis shown in Fig. [ija). We model the incident 
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optical axis 



Yz 



FIG. 2: Geometric relationship between the crystal coordi- 
nates (x, ct) and the optical axis coordinates {xo, cr) for both 
the reflected (top) and transmitted (bottom) wave. 



Equation ( 22 ) expresses the reflected wave from the tem- 
porally short incident field (21). Before attacking the 



integral, we must first relate the coordinates {x,t) to 
those along the new optical axis defined by the reflected 
wave. As depicted in Fig. [2] these comprise the coordi- 
nate transverse to the optical and the orthogonal 
time delay r defined at z = 0. Note that the time differ- 
ence is defined such that lines of constant t are parallel 
to the Xo-axis. Since the time t is measured with re- 
spect to lines of constant z, the reflected wave optical 
axis coordinates are related to {x,t) via x — Xo/ siii9 
and ct = CT + Xo/ tan 9 as shown in Fig. [2] 

To make further analytic progress, we can approximate 
this integral using the method of stationary phase. Com- 
putation of the stationary points is greatly simplified if 
we assume that |xo| ^ sin^ 6 and that the transverse size 
ax is sufficiently large. Specifically, we assume that 



wave by a Gaussian field both longitudinally and trans- 
versely that propagates along the optical axis defined by 
i sin 61 -I- xcosO: 
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If the incident pulse length ar is much shorter than the 
penetration length, we can approximate the temporal 
profile as a delta-function; taking the limit ct^ — J- and 
using the definitions ([5]), ( [T2| ), and (131, the incident field 
along the crystal surface z = in terms of the character- 
istic coordinates is 
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Upon inserting the initial condition (21) into the solu- 



tion for Ah (17), the integral over C can be trivially per- 



formed. To get the physical reflected field Eh requires 
the inverse Fourier transform as indicated by ([5|. The 
transform with respect to g is a Gaussian integral that 
can be taken analytically, so that 



Eh 
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Ji [vr sin 9{ct - x' cos 9) /A 
TT sin 6'(ct — x' cos 9) /A 
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(23) 



Physically, the first condition means that we consider 
distances behind the incident pulse cr that are much less 
than the Rayleigh range associated with the beam spot 
size projected onto the crystal surface, while the second 
condition is equivalent to requiring the beam angular di- 
vergence is much less than sin'^ 9; these conditions are 
typically well-satisfied for all but extreme grazing angles. 
Under these assumptions, the stationary point is 



X — ct cos ( 



sm 9 sm V sm 

and the reflected wave is given by 



P ikoXh „»yocfcor/2sin^ 9 '^1 [ttct/ (A sin 6>)] 
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From ( 25 ) , the temporal profile of the refiected wave 



oscillates over the time-scale (A/c) sm.9 according to the 
Bessel function Ji ; we graph the seguence of power peaks 
associated with \Ji{y)/y\ in Fig. ^a). Additionally, the 
field profile is laterally displaced along Xo as the time de- 
lay increases according to Xo = ct cot 9. We plot a spe- 
cific example of this behavior in Fig.|3][b) for the C(004) 
reflection at 56.85°, for which A w 22.8 /im and Aq ~ 1.5 
A. Fi gure|3|b) plots the reflected amplitude \Eh\ to more 
clearly identify the trailing pulses, which are displaced 
as predicted by Eq. (25). In the next section we apply a 



similar analysis to forward scattering. We will show that 
while the behavior of the transverse envelope is quite sim- 
ilar, the temporal profile of Bragg forward scattering is 
distinct, in that it depends importantly on the crystal 
thickness d. 



5 




FIG. 3: (a) The temporal power profile ~ \ Ji{y)/y^ of the Bragg reflected wave as a function of argument 7rcr/(AsinS). (b) 
Reflected field magnitude \Eh\ from the (004) Bragg reflection of diamond at an incidence angle 9 = 56.86° when the incident 
wave has ~ 10 fim. The wave is displaced transversely to the optical axis by an amount proportional to the time elapsed, 
with Xo = CT cot 6. 



B. Bragg forward scattering 

Having found the spatiotemporal dependence of the re- 
flected wave, we now turn to obtaining the transmitted 
wave Eq at the rear surface of the crystal. The solution 
can again be determined by Riemann's method; restrict- 
ing ourselves to d/ sinfl < ct < 3c?/ sin0, in the Appendix 
we show that the transmitted wave is given by 
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where the Riemann function that satisfies the relevant 



boundary conditions ( A12 ) and the adjoint equation ( 18 ) 

is m\ 
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Along the exit crystal surface CP we have ^ = C 
that the second integrand in ( 26 ) involves 
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The first term in ( 26 ) is merely the initial condition 
evolved along the characteristics, which is not the focus 
of this study, and we will henceforth drop it from our 
expressions. Of the remaining two terms, the one on the 
second line in ( 26 ) is generated by the reflected wave that 
is directly excited by the incident pulse, which can be 
shown to dominate the second term from the first line in 
the following way. Since Aq oc S{(' — x cos6 sm9/2), upon 
integration the second line scales as the product of 1 / A^ 
and the x-ray path length through the crystal d/ sm9. 
Using the expressions (17) and (21) for Ah, it is easily 



shown that the integral on the first line in ( 26 ) scales as 
l/A^ times a highly oscillatory function integrated over 
the time ct; for the times ct < d/ sin0 that we are con- 
sidering, the first line of ([26]) is therefore negligible with 
respect to that of the second line. 



Thus, the transmitted wave of interest is given by the second line in (26). An expression for the transmitted electric 
field envelope is obtained by proceeding as we did for the reflected wave; we use the definitions of the reduced fields 
(12)-([l3|, insert the function (28), and use the initially short incident pulse (21). The integration over (' is then 



trivial, and we find 



£o\p — 



smW 
An 



dx 



-iqx 



■ sin^ f /4o-^ g-i/co(a-Xo/2)(2C/ sin O+d/ sin 6 



-x' cos 9) 



dRo 



(29) 



C' = (2;'/2)cos« 



While this expression is difficult to interpret, we can make further progress by considering the field Eq in physical 
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space. We apply J dqe^'^^ to both sides of (291; again, the integral over q is a Gaussian one that can be taken 
analytically, leading to 
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Finally, the expression pO| ) should be written in term of the coordinates along the optical axis illustrated in Fig. [2j 
where Xo again labels the transverse position with respect to the transmission optical axis, while r is the relative 
time along the axis. From the Figure, we see that at the rear surface z — d we have x = {xq + dcos,0)/ sm0, and 
2C,/ sin0 = ct — d/ am0 — ct + x^j tan6'; these coordinates will be used in the final results. 

I 



While the expression ( 30 ) doesn't look any more appeal- 



ing than ( 29 ) , we can begin to make sense of it by consid- 
ering the limit of exact backscattering, i.e., — > 7r/2, in 
which case the integral over x is merely the Gaussian in- 
tegral associated with paraxial evolution. Furthermore, 
the transverse crystal coordinate equals its optical axis 
counterpart x = Xq^ while 2C, — ct. In the limit of exact 
backscattering ~ 7r/2, (30 1 simplifies to 
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In this form, the output field has a simple physical 
interpretation. The first factor shows that the ampli- 
tude scales as the product of the propagation distance 
and square of the coupling tt/A, since the incident wave 
must effectively be first scattered into the reflected wave 
and then back to the transmitted to produce the trail- 
ing pulse. Additionally, the amplitude is altered by both 
the phase difference and loss due to the crystal index 
of refraction and by the natural transverse spreading of 
the pulse: in the limit /cof^ 3> (d -|- ct), namely, that 
the Rayleigh range is much greater than the total prop- 
agation distance d -f cr, it yields the transverse envelope 
e~^°/'*'^^. Finally, the second line gives the transmission 
temporal profile or the "wake" of the incident pulse. 

We plot the temporal power profile ^ \Ji{v)/y\^ as a 
function of ttct/A for two crystal thicknesses in Fig.|4|^a). 
Doubling the crystal thickness from 5A to lOA decreases 
the time between maxima while increasing the peak 
power of each pulse by a factor of four. Note that the 
entire domain of the graph lies in the region of the first re- 
flection maximum plotted in Fig[3j^a); the characteristic 
time scale of the forward-scattered wave is much shorter 
than that of the reflected wave for d > A. Addition- 
ally, the temporal profile of forward scattering depends 
strongly on the crystal thickness d. 

Equation ( |31[ ) implies that the successive power max- 
ima of the transmitted field are given by the maxima of 
the function \ Ji{y)/y\, which correspond to the zeroes of 



the Bessel function J2{y)- We denote the positions of 
the \Jiiy)/y\ maxima by ^„ for integer n > 0, with the 
n — peak at — and the first trailing maximum 
of |Ji(y)/y| given by y = Setting the argument 

TT^CT{2d + ct)/JS? = Jl, we find 



C-(0 = 7r/2) 



d2 



7r2 
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(32) 



if the crystal thickness is much greater than the extinc- 
tion length, d 3> A. The time delay of the power maxima 
( [32| ) are inversely proportional to the crystal thickness d. 
This can be understood in a heuristic manner by con- 
sidering Fig. Qb). An initially sharp incident wave Eq, 
indicated in Fig. |4]jb) by the thick, red vertical arrow at 
C = 0, generates a reflected wave shown schematically as 
the blue horizontal lines directed along the ^ character- 
istics. As the reflected wave amplitude grows, it in turn 
couples to the transmitted wave (the red dotted lines), 
which increases its amplitude as it extracts energy from 
E^. In this way, energy alternates between the reflected 
and transmitted waves, as evidenced by the oscillatory 
profiles (251 and (31). The position of the first trans- 



mission maximum, for example, should be given by the 
integrated amplitude of the reflected wave just behind 
the incident pulse, which is proportional to the interac- 
tion volume given by the green shaded area in Fig. Qb) . 
If we denote the flrst transmission maximum as t^^^ , we 
have dr™^'^ = constant or r™^'^ ~ l jd. 

Equation (31) is rigorously valid but limited Xo — 
7r/2. We derive an approximate expression for (30) 



that is valid for arbitrary excluding grazing incidence 
by again using the method of stationary phase. First, 
we insert the optical axis coordinates shown in Fig. [2] 
which along the exit surface yield the replacements x — 
{xo + dcos0)/ sin 6* and 2C,/ sm0 = ct + Xo/tan6'. Next, 
we again assume that |xo| ^ sm^0, and that the trans- 
verse spot size is sufficiently large so that the Rayleigh 
range associated with the projected beam size is much 
shorter than the propagation length through the crystal, 
fcocr2sin^0 » {d/ sm0 + ct). Finally, we assume that 
the transverse size of the beam projected onto z is much 
smaller than the crystal thickness, so that jcTj; cos 6*1 <^ d. 
These conditions typically apply for all but small grazing 
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FIG. 4: (a) Time dependence of the Bragg forward diffraction intensity for a crystal thickness that is 5 and 10 extinction 
lengths. The position of first intensity maximum scales as 1/d while its peak power oc (d/ smOY . Note that in both cases there 
are several transmission maxima over a time interval for which Fig. |3|a) shows only one single reflection peak, (b) Illustration 
of the coupled forward and reflected modes in Bragg diffraction. The crystal surfaces are shown as dotted lines with unit slope. 
The initial pulse excites a reflected wave (blue arrows) that propagates along the C, characteristic, which in turn excites forward 
scattered waves (dotted red arrows) that move along ^. The first transmission maximum is determined by the coupling area 
shaded green, meaning that the time delay of the first maximum scales inversely with d. (c) Time dependence of Bragg forward 
diffraction intensity, plotted as a function of t/ sin 6* for the (004) reflection in diamond with d — 0.15 mm. The power maxima 
line up on the normalized t/ smd scale, with the peak scaling as (d/sinS)^. The central wavelength Ao ~ (1.79 A)sinS. 



incidence angles 9 <^1, and when satisfied result in a sta- 
tionary point that is identical to that of Bragg reflection 
(24): x'g = (xo — CT cot 9)/ sin 6*. Thus, the transmitted 
field amplitude is 



-^0 P — ^ 



h? sin 9 
X exp 



Jl 



[{xo - cr cot 6*)^] 
■K\JcT{2d/ sing + ct/ sin^ 9)/K 



(33) 



■KsJcT{2d/sai9 + cT/ sit? 9) /A 

From ( [33| , we see that the longitudinal time profile is 
nearly the same as that for exact backscattering ( [3l| ) 
with d — >■ d/ sui9 [this identification is exact for delays 
T <C {d/c)sm.9 if we assume that the Rayleigh range 
is sufficiently long]. Additionally, we recall that the ex- 
tinction length A is nearly independent of 0, so that the 
temporal envelope ^ Ji{y)/y can be obtained from that 
evaluated aX 9 — n/2 Eq. (31) by replacing t ^ t/ sai9. 
Thus, the maxima at arbitrary 9 can be written in terms 
of those at 6* = 7r/2 Eq. ((32| via 



T--(0) = C-(^/2)sin( 



(34) 



A2 



(35) 



with, when d ^ ^„A/7r, 

™^)=2.^ cd- 

We show in Fig. |4|^c) the time response of Bragg for- 
ward diffraction at three different incidence angles from 
the (004) Bragg reflection in diamond. The crystal thick- 
ness has been flxed at d = 0.15 mm, and we plot the 



power profiles as a function of t/ sm9. The position of 
the intensity maxima is invariant on the normalized time 
scale, while the power scales as [d/ sin 6')^. 



The temporal profile is modified as indicated in (34), 
while ( [33| ) also shows that the transverse envelope is 
translated from the optical axis at a:o = by the amount 
cr cot 0. This shift can be related to the distance along 
Xo that the reflected wave travels during the time t. We 
show this geometrically in Fig. [5]ja): much like how the 
delay arises because of the induced Bragg scattering in 
the crystal, the representative rays are now directed along 
both z and a;, so that as energy oscillates between waves 
the field is displaced in Xq. Here, the relevant shift can be 
found by comparing the times r and transverse coordi- 
nates Xo after the rays propagate along the characteristics 
some fixed distance. From Fig.[5ja), we see that while the 
incident ray propagates a distance D, the refiected wave 
is displaced in time by an amount ct = D — D cos 29 
and transversely hy Xo = D sin 29. Taking the ratio, we 
find that corresponding to a time delay r, we have the 
displacement Xo — cr sin 26'/(l — cos 20) = cr cot 6*. 

We show the spatiotemporal profile of a Bragg forward 
scattered pulse from the C(004) crystal aX 9 — 56.86° in 
Fig. [5|b). We assume a temporally short input pulse 
with rms transverse width cr^. = 10 /im. Again, we plot 
the magnitude |i?o| to more clearly show all the trailing 
pulses, and we scale the amplitude so that the first de- 
layed maxima has unit magnitude. The temporal pro- 
file closely mirrors that shown of the same crystal in 
Fig. |4]jc), while each subsequent pulse is more distantly 
displaced from the optical = 0. The cyan line is 

drawn along the theoretically predicted line Xo = ct cot 0, 
which closely predicts the transverse displacement. In 
this case, the first trailing pulse is displaced from the op- 
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FIG. 5: (a) Illustration of the coupled forward and reflected modes in the crystal. As explained in the text, a forward 
scattered ray is displaced Aa; = D sin 25 laterally while delayed by the time r = D{1 — cos26)/c; eliminating D, we find that 
Aa; = CT cot d. (b) Spatiotemporal dependence of the electric field magnitude in forward Bragg diffraction. We chose the 
(004) reflection from a 0.1-mm-thick diamond crystal at an incidence angle 9 = 56.86° (Aq = 1.5 A), while the incident pulse 
has an rms width ~ 10 fj,m. The temporal dependence can be compared with that plotted in Fig. |4|^c), while the lateral 
displacement Axo = cr cot 6 is emphasized by the cyan arrow. 



tical axis by a small amount ^ 5 fim, while the second 
maxima is shifted by an amount ^ 13 /im, which is of 
the order of the rms width. 



II. MONOCHROMATIC POWER FOR FEL 
SELF-SEEDING 

We have seen that Bragg forward scattering gives rise 
to a sequence of delayed power maxima from temporally 
short incident radiation. Reference ^ proposed using 
the first trailing maximum to seed an FEL at hard x-ray 
wavelengths: an initially short SASE pulse generates the 
radiation "wake" depicted in Fig. [5jb) that is then used 
to coherently seed the FEL interaction in downstream 
undulators. In the preceding discussion, we computed 
the spatiotemporal field profile generated by an initially 
short and coherent incident pulse; since SASE is tempo- 
rally incoherent (chaotic light), however, some additional 
considerations are necessary to determine the relevant ra- 
diation seed power. 

The longitudinal structure of SASE can be well mod- 
eled as a sum of Gaussian modes that have random tem- 
poral positions and phases. For M longitudinal modes, 
we approximate the SASE field as 



M 

E 



(36) 



where U{t) is the total envelope of the field determined 
by the electron beam current and FEL gain, tj are a 
random set of times, and Cj are a set of random complex 



amplitudes such that M{\tjf') /ar is proportional to the 
ensemble-averaged SASE energy. The temporal width a-r 
is dictated by the FEL physics while the number of modes 
M is determined by ar and the characteristic length of 
the envelope U{t) that we denote Lpuisc- Assuming that 
CTr is sufficiently short, the incident pulse (36) generates 



a trailing electromagnetic field that is merely a sum of 
M wakes discussed previously, with each one beginning 
at the time tj and having the relative complex amplitude 
given by Ej. 

First, we consider the case when the duration of the 
incident SASE length Lpuisc is much shorter than the 
characteristic time-scale of the trailing wake, which is the 
typical situation for self-seeding of few-fs pulses. In this 
case, the M Bragg forward-diffracted beams generated 
by (36) have the same temporal shape, but differ with 



random phases and amplitudes. The ensemble-averaged 
total power in the trailing pulse scales as M times the 
average power in an individual wake. If we consider the 
power maxima of ( 33 ) , the ratio of the ensemble-averaged 



power in the first trailing seed to that in the SASE is 



pmax 
seed 



= A-KdiM 



pulse 



TT^d Ji( 

A2 sine 



(37) 



where again is the position of the first maxima of 
\Ji{y)ly\ with y > 0. The temporal width of each mode 
and number of modes is dictated by the FEL process, 
and as such will depend on the electron beam quality 
and profile. In order to get an approximate expression, 
we use the analytic results that are available for a long 
electron beam with zero initial energy spread [31 [27] . In 
this case, the rms width of a SASE mode is related to the 
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FEL coherence length Lcoh by Ut = _Lcoh/v 37r while the 
number of longitudinal modes for a pulse of total length 
-^^puise is M = Lpuisc/-^coh; inserting this into ^\ yields 



pmax 

seed 

-fpulsc 



4 



pulse 



n 2 



A2 sin 9 



(38) 



Finally, the SASE coherence length generated by an elec- 
tron beam of vanishing energy spread is 



coh 




(39) 



where Nq is the number of gain lengths in the upstream 
undulator. Combining (39) and (37 1, we have 



Pr 



pulse 




0.0058 



(40) 



Due to the various approximations made, we only ex- 



pect (40) to be correct to within a factor of two or so. 



Nevertheless, this yields a useful estimate and has the ap- 
propriate scaling when Lpuisc is much less than the time 
scale of variation of the forward Bragg diffraction signal; 
if we use (34)-(|35|) to eliminate the crystal parameters in 



favor of the delay, ( 40 ) becomes 




pmax 
seed 



pulse 



-^coh-^pulsc 

(cri)2 ■ 



(41) 



For example, for the approximate parameters proposed in 
[28] for use at the Linac Coherent Light Source (LCLS), 
the low charge operation has ipuigg ~ 1 /^^m and after the 
suggested interaction length we find that Lcoh ~ 0.07 /im; 
for a 0.1-mm diamond crystal at an angle 9 = 56.86°, 
the C(004) planes have n'^d/iA'^ sine) « 2.4 /im , and 
(40) implies that a 1-GW SASE pulse would produce 



a seed whose ensemble-averaged peak power is 2.3 MW; 
the same result obtains from ( 41 ) for the associated delay 
Ti « 18 fs. These estimates compare quite favorably with 
the quoted result of about 2.5 MW [21] . 



To demonstrate the scaling predicted in ( 40 ) holds in 
general (assuming that Lpuise 'C A^/ir^d), we used a sim- 
ple ID FEL code to generate SASE output for several 
different pulse lengths Lpuise and over two different num- 
ber of gain lengths to compute the seeding wake after the 
monochromator for LCLS-type parameters cited above. 
We show the results in Fig. [6j where we scale the seed 
power with the ratio (40), so if our results are exact all 
the lines should overlap with the theory line. We note 
that the seed power is correctly predicted within 10% or 
so, while the delay of the seed increases over that pre- 
dicted by the theory by an amount approximately equal 
to Lpuise- For short pulses, this shift is a minor correc- 
tion. 

The power statistics of the seeding wake are identical 
to that of a single SASE intensity mode/spike, meaning 



6 fs, Ng = 8.8 
6 fs, Ng = 13.5 . 
3 fs, Ng ■■ 
3 fs, Ng = 13.5 — 
2 fs, Ng = 13.5 — 
Theory 




FIG. 6: Ensemble-averaged power of the monochromatic seed 
obtained from simulation, scaled to the theoretical predicted 
amplitude ( 40 1 . We have varied both the electron beam tem- 



poral duration and the undulator length to produce the Lpuiac 
and Ng listed in the key. We have used a fiat-top electron 
beam profile to more accurately determine Lpuisc; more real- 
istic distributions can be expected to yield power deviations 
from the theory by a factor of two. The position of the first 
maximum linearly increases as the pulse duration increases, 
and therefore varies by a few femtoseconds. 



that the effective seed power fluctuates 100%. In the 
frequency domain, this reflects the fact that the power 
fluctuations of chaotic light within a bandwidth much 
less than c/ctt approach unity. 

Finally, we wish to make a few statements in the ap- 
plication of Bragg forward scattering to the seeding of 
relatively long pulses. By long pulses, we mean those for 
which Lpuiso ^ A^ sin6'/7r^d; for the LCLS parameters 
used herein, this applies for pulses longer than about 10 
fs. In this case, the temporal modes in (36) are spread 



over a time that is longer than that associated with Bragg 
forward scattering. Thus, we expect that the effective 
number of modes contributing to the seed amplitude to 
be M < (A^ sin9/'K^d)/ar- While the peak power in this 
trailing monochromatic seed is therefore of order that 
in the short pulse case, one must confront the fact that 
as the delay t is increased to accommodate the longer 
pulse, the wake is displaced transversely according to 
Xo = CT cot 9. Since a typical FEL-produced radiation 
beam has ^ 20 /im and cot 9 ^ 1 to allow variations 
in the central wavelength, this limits the applicability of 
this self-seeding scheme to Lpuiso/c ^ 50 fs. 



III. CONCLUSIONS AND FUTURE 
DIRECTIONS 



We have used the dynamical theory of x-ray diffraction 
to calculate approximate, analytic expressions for both 
the reflected and forward scattered electric field ampli- 
tudes resulting from Bragg diffraction of an initially short 
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x-ray pulse incident on a crystal. Both field profiles are 
characterized by a temporal envelope that oscillates in 
time ^ Jiiy)/y, where y oc t/ sin 9 for the reflected wave 
while y cx \J At j sin 9 for the forward diffracted field. 
The delayed output is also displaced transversely by an 
amount Xo = ct cot 6*, which can be associated with the 
coupled wave interaction in the crystal. Finally, we used 
the developed theory to analyze in detail the dynamics 
of the "wake monochromator" as applied to self-seeding 
for hard x-ray free-electron lasers. We found simple rela- 
tionship for the induced delay and output power in terms 
of the crystal and SASE parameters. 

In the future we plan to extend this treatment to asym- 
metrically cut crystals in both Bragg and Laue scattering 
geometries, which can be used to analyze the field output 
generated by temporally short x-ray pulses from a wide 
variety of x-ray optical elements. 
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Appendix A: Riemann's method and application to 
Bragg scattering 

Any linear second-order differential operator £ that is 
also hyperbolic can be written as acting on the function 
A via (see, e.g., ^5] ) 



«g^-4 + 6— ^ + cA (Al) 



where (C',C') ^-re the characteristic coordinates. We ad- 
ditionally consider the adjoint £^ of the operator £, 
which is defined through the inner product relation- 
ship J RC{A) = / C''{R)A. This definition implies that 
RC{A) — £^{R)A can be written as the divergence of a 
vector function 



RC(A)~L\R)A = V.P='^ + '^ 



(A2) 



with P vanishing at the endpoints. Thus, is com- 
pletely specified by both the equation ( |A2p and by a set 
of boundary conditions. Ignoring for the moment these 
boundary conditions, one can easily show that the differ- 
ential identity ( A2 1 is satisfied if the action of adjoint 
is given by [29] 



d_ 



-{bR)+cR^O, (A3) 



and the components of P 



dA 



Pc = 



R 



R 



dA 



OR 
dC 
dR 



are 



A- 



-A 



+ aRA 



bRA. 



(A4) 



(A5) 



Riemann's met hod is obtained by integrating the differ- 
ential identity ( A2 1 over an arbitrary closed region F in 
the (C'jC) plane. Using Green's identity on the right- 
hand side, we find 



J da [RCiA) - AC^R)] = jdi 7 • P, 



(A6) 



where the integration proceeds along the boundary 7 and 
7 is its outward-facing normal. Now, if we assume that 
A \& a, solution to C{A) — 0, and R solves the adjoint 
equation C^{R) — 0, then fd£ 7-P vanishes. By inserting 
P, we have an expression that can be used to find the 
solution to C{A) = (assuming it exists) in terms of its 
values along the boundary and the Riemann function R. 

For the problem of Bragg scattering, a = 6 = 0, so 
that C and satisfy the same partial differential equa- 
tion, and C is formally self-adjoint. Now, we can find the 
solution by determining an appropriate curve 7 and Rie- 
mann function R which, due to the differing boundary 
conditions, will depend on whether we consider forward 
or backward scattering. In the latter case, the triangle 
ABQ is a convenient choice to obtain the reflected wave 
Ah at point Q, which is chosen such that < mean- 
ing that the solution is given only over the time interval 
< t < 2d/sm9. If we further assume that the fields 



Ao, Ah vanish along the line AB, we can integrate (A6) 
by parts to write it as 





Q 


RAh\Q- 


/' 




B 








Ah 





dC Ah 
dR 



dR 



(A7) 



R 



dAh 
dC 



Here, the Riemann function R — R{C,£,;C\£,') satisfies 
the adjoint equation associated with the (formally self- 
adjoint) system (14): 



92 



d^'dC 



7^(C,e;C',r) 



A2 



R{c^; ea- 



rn 



If R additionally satisfies dR/d^' = along AQ (the 
line C' = C) and dR/dC,' = along liQ (when ^ 



then (A7| gives the solution for the reflected wave in 
terms of the incident wave Aq oc dAh/dC^' . The Rie- 
mann function that satisfies these requirements and also 
i?(g) = i?(C,f;C',e')-lis[22l 



2vrV(C-C')(C-e')/A 



(A9) 
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For the Riemann function Rh only the integration along 
AQ contributes from ( A7 1 , and the field is therefore given 
by 



dC 't^MC, ORkiC C, O- (Aio) 

z sm u 



To determine the equation for the transmitted wave at 
point P on the rear surface of the crystal CP, we choose 
the contour 7 to be the parallelogram ACPQ, which is 
restricted to d/ ainO < t < 3d/ sin 6* to exclude the effects 
of surface refiections. Again assuming that the fields ^0 
and Ah vanish along the line AB, (A6I can be partially 
integrated to yield 



^ , , OR dAo 

Q 



j^dAo OR 



(All) 



Along the front surface the incident wave Aq is prescribed 
while the reflected wave Ah i s given by (AlO). On the 



other hand, the reflected wave Ah oc dAa/dSJ vanishes 
along the rear surface CP. Thus, to reduce (All) to a 



closed form solution, we require that the Riemann func- 
tion i?o solving (A8) also satisfles the auxilliary condi- 
tions 



dRo 



PQ 



dRo 



0, Ro{P) = 1. (A12) 



CP 



The Riemann function satisfying the boundary condi- 
tions (A12| and the adjoint equation (A8) associated 



with the transmitted wave is 



27rV(C-C')(e-e)/A' 



(A13) 



while rearranging (All I implies that the solution 



A|p = i?o-4oU- /<i?o(C,C;C',C')^-4 

' 2 sm d 



dC Ao{C,C)g^Ro{C,^;e,C) 
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